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Abstract 



We map the general relativistic two-body problem onto that of a test parti- 
cle moving in an effective external metric. This effective-one-body approach 
defines, in a non-perturbative manner, the late dynamical evolution of a coa- 
lescing binary system of compact objects. The transition from the adiabatic 
inspiral, driven by gravitational radiation damping, to an unstable plunge, 
induced by strong spacetime curvature, is predicted to occur for orbits more 
tightly bound than the innermost stable circular orbit in a Schwarzschild 
metric of mass M = mi + m<i- The binding energy, angular momentum and 
orbital frequency of the innermost stable circular orbit for the time-symmetric 
two-body problem are determined as a function of the mass ratio. 



I. INTRODUCTION 



Binary systems made of compact objects (neutron stars or black holes), and driven toward 
coalescence by gravitational radiation damping, are among the most promising candidate 
sources for interferometric gravitational- wave detectors such as LIGO and VIRGO. It is 
therefore important to study the late dynamical evolution of a coalescing binary system of 
compact objects and, in particular, to estimate when occurs the transition from an adiabatic 
inspiral, driven by gravitational radiation damping, to an unstable plunge, induced by strong 
spacetime curvature. The global structure of the gravitational wave signal emitted by a 
coalescing binary depends sensitively on the location of the transition from inspiral to plunge. 
For instance, in the case of a system of two equal-mass neutron stars, if this transition occurs 
for relatively loosely bound orbits, the inspiral phase will evolve into a plunge phase before 
tidal disruption takes place. On the other hand, if the transition occurs for tightly bound 
orbits, tidal effects will dominate the late dynamical evolution. 

In this paper we introduce a novel approach to the general relativistic two-body problem. 
The basic idea is to map (by a canonical transformation) the two-body problem onto an 
effective one-body problem, i.e. the motion of a test particle in some effective external metric. 
When turning off radiation damping, the effective metric will be a static and spherically 
symmetric deformation of the Schwarzschild metric. [The deformation parameter is the 
symmetric mass ratio v = m 1 m 2 / '{mi + m 2 ) 2 .] Solving exactly the effective problem of 
a test particle in this deformed Schwarzschild metric amounts to introducing a particular 
non-perturbative method for re-summing the post-Newtonian expansion of the equations of 
motion. 

Our effective one-body approach is inspired by (though different from) an approach 
to electromagnetically interacting quantum two-body problems developed in the works of 
Brezin, Itzykson and Zinn- Justin |l| (see also ||) and of Todorov and coworkers [[|, |4[]. 
Ref. has shown that an approximate summation (corresponding to the eikonal approx- 
imation) of the "crossed-ladder" Feynman diagrams for the quantum scattering of two 
charged particles led to a "relativistic Balmer formula" for the squared mass of bound states 
which correctly included recoil effects (i.e. effects linked to the finite symmetric mass ratio 
v = m 1 m 2 /(m 1 + m2) 2 ). However, the eikonal approximation does not capture some of the 
centrifugal barrier shifts which have to be added by hand through a shift n —>■ n — ej of 
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the principal quantum number JT]], ||. The approach of Ref. J3| is more systematic, being 
based on a (Lippmann-Schwinger-type) quasi-potential equation whose solution is fitted to 
the Feynman expansion of the (on-shell) scattering amplitudes {p'ip'2 |5'|piP2}- However, 
several arbitrary choices have to be made to define the (off-shell) quasi-potential equation 
and the nice form of the relativistic Balmer formula proposed in Ref. |l| is recovered only 
at the end, after two seemingly accidental simplifications: (i) the ratio of some compli- 
cated energy- dependent quantities simplifies || , and (ii) the second-order contribution to the 
quasi-potential contributes only to third order. We note also that the extension of Todorov's 
quasi-potential approach (initially developed for quantum two-body electrodynamics) to the 
gravitational two-body problem [f§] leads to much more complicated expressions than the 
approach developed here. 

Before entering the technical details of the effective one-body approach, let us outline 
the main features of our work. We use as input the explicit, post-Newtonian (PN) expanded 
classical equations of motion of a gravitationally interacting system of two compact objects. 
In harmonic coordinates (which are convenient to start with because they are standardly 
used for computing the generation of gravitational radiation), these equations of motion 
are explicitly known up to the 2.5PN level ((w/c) 5 -accuracy) ||, 0. They have the form 
(a, b =1,2) 

a a = Af N {z b , v h ) + Ar c (z b , v b ) + 0(c~ 6 ) , (1.1) 

where *4. 2PN = Aq + c~ 2 A2 + c~ 4 A4 denotes the time-symmetric part of the equations 
of motion, and A reac = c -5 A 5 their time-antisymmetric part. Here, z a , v a , a a , denote 
the positions, velocities and accelerations, in harmonic coordinates, of the two bodies. [In 
this work we consider only non-spinning objects.] Throughout this paper, we shall use the 
following notation for the quantities related to the masses m x and m2 of the two bodies: 

B/f _ . _m l m 2 _ \i _ m x m 2 

M = m l + m 2 , » = "=M= {mi+m2) 2 - 

Note that the "symmetric mass ratio" v varies between (test mass limit) and \ (equal mass 
case). 

We first focus on the time-symmetric, 2PN dynamics defined by J\^ N (z b , v b ). After 
going to the center of mass frame (uniquely defined by the Poincare symmetries of the 
2PN dynamics), and after a suitable coordinate transformation (from harmonic coordinates 
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to ADM coordinates z a — > q a ), the dynamics of the relative coordinates q = q l — q 2 is 
defined by a 2PN Hamiltonian H(q,p). Starting from H(q,p), we shall uniquely introduce 
a 2PN-accurate static and spherically symmetric "effective metric" 

ds 2 cS = —A(R cS ) c 2 dt 2 ca + dR *z + ^(^eff + sin 2 9 cS d^) , (1-3) 

where 

such that the "linearized" effective metric (defined by a\ and d\) is the linearized 
Schwarzschild metric defined by the total mass M = mi + m 2 , and such that the effective 
Hamiltonian H e s(q eS ,p eS ) defined by the geodesic action — J /icds e g, where \i = m\m2jM 
is the reduced mass, can be mapped onto the relative-motion 2PN Hamiltonian H(q,p) by 
the combination of a canonical transformation (q e s,p e g) — > (Q,p) and of an energy trans- 
formation H = f(H e ff), corresponding to an energy- dependent "canonical" rescaling of the 
time coordinate dt e g = dt (dH / dH e s) . 

The effective metric so constructed is a deformation of the Schwarzschild metric, with 
deformation parameter the symmetric mass ratio v = fi/M. Considering this deformed 
Schwarzschild metric as an exact external metric then defines (in the effective coordinates) 
a //-deformed Schwarzschild-like dynamics, which can be mapped back onto the original 
coordinates q a or z a . Our construction can be seen as a non-perturbative way of re-summing 
the post-Newtonian expansion in the relativistic regime where GMj {c 2 \q 1 — q 2 \) becomes of 
order unity. In particular, our construction defines a specific //-deformed innermost stable 
circular orbit (ISCO). Superposing the gravitational reaction force A reac onto the "exact" 
deformed-Schwarzschild dynamics (defined by mapping back the effective problem onto the 
real one) finally defines, in a non-perturbative manner, a dynamical system which is a good 
candidate for describing the late stages of evolution of a coalescing compact binary. 



II. SECOND POST-NEWTONIAN DYNAMICS OF THE RELATIVE MOTION 

OF A TWO-BODY SYSTEM 



Let us recall some of the basic properties of the dynamics defined by neglecting the 
time-odd reaction force in the D amour- Deruelle equations of motion ( |1.1|) . The 2PN (i.e. 
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(v /c) 4 -accurate) truncation of these equations of motion defines a time-symmetric dynamics 
which is derivable from a generalized Lagrangian L(z±, z 2 , Vi, v 2 , Oi, a 2 ) ||, J7| (a function 
of the harmonic positions, z±, z 2 , velocities v±, v 2 and accelerations a±, a 2 ). The generalized 
Lagrangian L(z\, z 2 , vi, v 2 , a\, a 2 ) is (approximately) invariant under the Poincare group 
0. This invariance leads (via Noether's theorem) to the explicit construction of the usual 
ten relativistic conserved quantities for a dynamical system: energy S, linear momentum "P, 
angular momentum JT, and center-of-mass constant K, — Q — 'P t. Because of the freedom 
to perform a Poincare transformation (in harmonic coordinates), we can go to the (2PN) 
center-of-mass frame, defined such as 



T = K = Q = 0. (2.1) 

Refs. ||10|| , |TT| explicitly constructed the coordinate transformation between the harmonic 
(or De Donder) coordinates, say used in the Damour-Deruelle equations of motion, and 



the coordinates, say g M , introduced by Arnowitt, Deser and Misner [12] in the framework of 
their canonical approach to the dynamics of the gravitational field. The Lagrangian giving 
the 2PN motion in ADM coordinates has the advantage of being an ordinary Lagrangian 
L(q x , q 2 ,<li,q 2 ) (depending only on positions and velocities), which is equivalent to an ordi- 



nary Hamiltonian H(q ± , q 2 ,Pi,p 2 ) ]I3" |, ||14j| . The explicit expression of the 2PN Hamiltonian 
in ADM coordinates, H(q ly q 2 ,p 1 ,p 2 ), has been derived in Ref. |L1| by applying a contact 
transformation 

q a (t) = z a (t)-5*z a (z,v), (2.2) 

to the generalized Lagrangian L(z a , v a , a a ). The shift 5* z a is of order 0(c~ 4 ) and is defined 
in equation (35) of [jlO[ or equations (2.4) of The contact transformation ( |2.2| ) removes 



the acceleration dependence of the harmonic-coordinates Lagrangian L harm (z, v, a) and trans- 
forms it into the ADM-coordinates ordinary Lagrangian L ADM (q, q). A further Legendre- 
transform turns L ADM (g 1; q 2 , q ± , q 2 ) into the needed 2PN Hamiltonian H(q 1) q 2l p 1) p 2 ) in 
ADM coordinates. The explicit expression of this Hamiltonian is given in equation (2.5) of 
Ref. ]TTJ. It has also been shown in Ref. |IJ that the Hamiltonian H(q l ,q 2 ,p l ,p 2 ) can be 



directly derived in ADM coordinates from the (not fully explicit) iV-body results of Ref. [13 



by computing a certain integral entering the two-body interaction potential. [For further 
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references on the general relativistic problem of motion, see the review [15|]; for recent work 
on the gravitational Hamiltonian see |]16| , ]T7| , [[TJ|, 



The ADM expression of the total Noether linear momentum "P associated to the trans- 
lational invariance of L(z,v,a) is simply "P = p x + p 2 . Therefore it is easily checked that, 
in the center-of-mass frame (|2.1| ), the relative motion is obtained by substituting in the 
two-body Hamiltonian H(q 1 ,q 2 ,p 1 ,p 2 ), 

Pl ^P, P 2 ^-P, (2.3) 

where P = dS/dQ is the canonical momentum associated with the relative ADM position 
vector Q = q 1 — q 2 . [For clarity, we modify the notation of Ref. [Tl[] by using q l7 q 2 , Q 
and q for the ADM position coordinates which are denoted r 1 , r 2 , R and r , respectively, 
in Ref. iTl.1 



Our technical starting point in this work will be the reduced center-of-mass 2PN Hamil- 
tonian (in reduced ADM coordinates). We introduce the following reduced variables (all 
defined in ADM coordinates, and in the center-of-mass frame): 

= Q_ = <?i ~ 9 2 = P 

~ GM ~ GM ' P ~ fi ' 

t ~ # NR H R -Mc 2 , . 

In the last equation, the superscript "NR" means "non-relativistic" (i.e. after subtraction of 
the appropriate rest-mass contribution), while "R" means "relativistic" (i.e. including the 



appropriate rest-mass contribution). From equation (3.1) of JTT] the reduced 2PN relative- 
motion Hamiltonian (without the rest-mass contribution) reads 

H(q,p) = H (q,p) + -4 H 2 (q,p) + \ H 4 (q,p) , (2.5) 



where 



& c 4 



H (q,p) = l-p 2 --, (2.6a) 
L q 

H 2 (q, p) = ~ (1 - 3u) p A - 1 [(3 + v) p 2 + u(n ■ p) 2 } + ^ , (2.6b) 

H A (q,p) = ^-(l-5u + 5u 2 )p e 
16 

+ — [(5 - 20z/ - 3z/ 2 ) p 4 - 2u 2 p 2 (n ■ p) 2 - 3u 2 (n ■ p) 4 ] 
8q 

+ ^ [(5 + 8u) p 2 + 3u(n- p) 2 } - ± (1 + 3v) , (2.6c) 
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in which q = \q\ = (q 2 ) 1 ^ 2 and n = q/q. When convenient, we shall also use the notation r 
for the reduced radial separation q (and R for the unreduced one Q). [As in Eqs. ( |2.8| )- (|2.12| ) 
below. 1 



The relative-motion Hamiltonian (|2.5|) is invariant under time translations and space 



rotations. The associated conserved quantities are the reduced center-of-mass (cm.) energy 
and angular momentum of the binary system: 

H(q lP ) = S* R = ^ } qx P = 3 = ^. (2.7) 

A convenient way of solving the 2PN relative-motion dynamics is to use the Hamilton- Jacobi 
approach. The motion in the plane of the relative trajectory is obtained, in polar coordinates 

q x = rcosip, q y = rsimp, q z = , (2.8) 

by separating the time and angular coordinates in the (planar) reduced action 

S = = -£* R t + j cp + S r (r, £* R , j) . (2.9) 

The time- independent Hamilton- Jacobi equation H NR (q,p) = S^ R with p = dS/dq can be 
(iteratively) solved with respect to (dS r /dr) 2 with a result of the form 

S r (r,S^ R ,j) = Jdr y/K(r,&™,j). (2.10) 

The radial "effective potential" lZ(r, ^ NR ,j) is a fifth-order polynomial inl/r = 1/q which 
is explicitly written down in equations (3.4) of JTT]. In this section, we shall only need the 



corresponding (integrated) radial action variable 

. Ir 2 



fiGM 2tt 



/Tmai I 

/ dr^/n(r,£^,j). (2.11) 



The function i r (E^ K , j) has been computed, at the 2PN accuracy, in Ref. |TT] (see equation 
(3.10) there). To clarify some issues connected with the fact that the natural scalings in the 
"effective one-body problem" (to be considered below) differ from those in the present, real 
two-body problem, let us quote the expression of the unsealed radial action variable 

I R = ai r = — / dR — , (2.12) 

2tt JR mia dR 
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in terms of the unsealed variables £ NR = /x £ NR and J = aj. Here R = Q = GMr = GMq, 
and we introduced the shorthand notation 



a = ji GM = Gmi m 2 
for the gravitational two-body coupling constant. We have 



(2.13) 



a /i 



1/2 



V-2£ NR 
- J + 



Q- 



C 2 J 



15 

T 



, 15 



" 4 

3z/ 



5 



NR 



gNR 
/IC 2 



+ 



35 15 

1 v 

32 16 



35 
T 



5 

2" 



32 



a 



z/ 



NR 



/iC z 



(2.14) 



Equation ( |2.14j ) can also be solved with respect to £ NR = £ R — Mc 2 with the (2PN-accurate) 
result (see equation (3.13) of Ref. [|TT| ) 



£K { M,J) = Mc>-\^ 



1 + 



Q- 



6 



1 15 -v 



or i 5 7 

I 7 ^ 



2i/ 



+ 



27 



c 4 V2 AfJ 3 AT 2 J 



HJ AM 2 

3 35 - Au 1 145 - Ihv + z/ 2 
~ 2 AT 3 J + 8 



(2.15) 



where M denotes the Delaunay action variable M = Ir + J . The notation is chosen so as 
to evoke the one often used in the quantum Coulomb problem. Indeed, the classical action 
variables Ir and J7, or their combinations M = Ir + J and J7, are adiabatic invariants 
which, according to the Bohr-Sommerfeld rules, become (approximately) quantized in units 
of h for the corresponding quantum bound states. More precisely N / h becomes the "principal 
quantum number" and J jft the total angular- momentum quantum number. The fact that 
the Newtonian- level non-relativistic energy £ NR — — | jia 2 /N 2 + O(o~ 2 ) depends only on 
the combination M = Ir + J is the famous special degeneracy of the Coulomb problem. 
Note that 1PN (and 2PN) effects lift this degeneracy by bringing an extra dependence on 
J . There remains, however, the degeneracy associated with the spherical symmetry of the 
problem, which implies that the energy does not depend on the "magnetic quantum number" , 
i.e. on Ai = J z , but only on the magnitude of the angular momentum vector J = a/ J' 2 . 
Though we shall only be interested in the classical gravitational two-body problem, it is 
conceptually useful to think in terms of the associated quantum problem. From this point of 
view, the formula ( |2.15 ) describes, when Af/h and J /h take (non zero) integer values, all the 
quantum energy levels as a function of the parameters M = m x +m 2 , /i = mi m 2 /(m 1 +m 2 ), 
a = Gm\m 2 and v = fi/M. It is to be noted that the function £ R (Af, J) describing the 
energy levels is a coordinate- invariant object. 
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III. SECOND POST-NEWTONIAN ENERGY LEVELS OF THE EFFECTIVE 

ONE-BODY PROBLEM 



The "energy levels" (|2.15|) summarize, at the 2PN accuracy, the dynamics obtained by 



eliminating the field variables g^x) in the total action of a gravitationally interacting binary 
system 

Stot[£i,4\sv] = - J m 1 cds 1 - J m 2 cds 2 + Sfi e id[5V(z)] , (3.1) 



where ds\ = y—g^jjz\) dz^ dz\ and where S^sidig^x)] is the (gauge-fixed) Einstein-Hilbert 
action for the gravitational field. Let S TeB j[zi, z%\ be the Fokker-type action obtained by 
(formally) integrating out g^ v {x) in (|3.1|). [See, e.g., [|It]] for more details on Fokker-type 



actions. As we work here only at the 2PN level, and take advantage of the explicit results 
of Refs. @] , 0, we do not need to enter the subtleties of the elimination of the field degrees 



of freedom, which are probably best treated within the ADM approach. See 0], ||14|| .1 

The basic idea of the present work is to, somehow, associate to the "real" two-body 
dynamics S Tea \[zi , z%] some "effective" one-body dynamics in an external spacetime, as de- 
scribed by the action 

S eS [z%] = - m cds , (3.2) 



where dso = J— g^^o) dz^dz^, with some spherically symmetric static effective metric 



ds 2 eS = gf u (x x eS ) dx^ d< ff = -A(R cS ) c 2 dt 2 cS + B(R cS ) dR 2 eS 

+ C(R cS ) R 2 cS (d6 2 cfi + sin 2 9 cS d<p 2 cS ) . (3.3) 

To ease the notation we shall, henceforth in this section, suppress the subscript "eff" on the 
coordinates used in the effective problem. [Later in this paper we shall explicitly relate the 
coordinates Zq of the effective particle to the coordinates z±, z% of the two real particles.] 
The metric functions A(R), B(R), C(R) will be constructed in the form of an expansion in 
1/R: 

A{R) = 1 + ^R + ^k + ^k + '" ' 

B ^ = 1 + )R + Jh + -- ^ 
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Beware that the variable R in Eqs. (|3.4j) denotes (in this section) the effective radial coordi- 
nate, which differs from the real ADM separation Q = -Radm = GMr used in the previous 
section (e.g. in the definition of Ir). We indicate in Eq. ( |3.4| ) the terms that we shall need 
at the 2PN level. The third function C(R) entering the effective metric will be either fixed 
to Cs(R) = 1 (in "Schwarzschild" coordinates), or to satisfy Cj(R) = B(R) (in "Isotropic" 
coordinates) . 

There are two mass parameters entering the effective problem: (i) the mass m of the 
effective particle, and (ii) some mass parameter Mo used to scale the coefficients a^, bi entering 
the effective metric. For instance, we can define Mq by conventionally setting 

ai = -2GM . (3.5) 

By analogy to Eq. (|2.15| ), we can summarize, in a coordinate-invariant manner, the 
dynamics of the effective one-body problem ( |3.2| )-( |3~4| ) by considering the "energy levels" of 
the bound states of the particle m in the metric g e ^ u \ 

£R = m c 2 + £ NR = F{N , Jo; mo, a u h) . (3.6) 

Here, the relativistic effective energy £^ and the effective action variables Mq, Jo are unam- 



biguously defined by the action (|3.2|) . Namely, we can separate the effective Hamilton- Jacobi 
equation 

aSc£ aSeff 2 2 _ , . 

9 eS dx » QxV + m oC -U, (3.0 
by writing (considering, for simplicity, only motions in the equatorial plane 9 = | ) 

S cS = -£ t + J <p + S° R (R, So, Jo) ■ (3.8) 

To abbreviate the notation we suppress the superscript "i?" on the relativistic effective 
energy So- Inserting Eq. ( |3.8| ) in Eq. ( |3.7| ) yields 

1 Sq 1 f dS R \ J 2 



+ -57^hj# +7^ + m o^ = 0, (3.9) 



A(R) c 2 B(R) \dR J C(R)R 2 

and therefore 



S R (R,£ ,J ) = j dR^Ko(R,S ,J ), (3.10) 
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where 



TZq(R, £q, Jo 



B(R) gp 
A(R) c 2 



B(R) " 2 



m c 



Jo 



C(R) R 2 



(3.11) 



The effective radial action variable I R is then defined as 



Ir(£o, Jo) = 7T V^o{R, £o, Jo) 

Zn 



(3.12) 



Rn 



while the effective "principal" action variable Ao is defined as the combination J\f = I r + Jq. 

To obtain the effective "energy levels" £o = J- (No, Jo) one needs to compute the definite 
radial integral (|3.12 ). Ref. |11| (extending some classic work of Sommerfeld, used in the old 



quantum theory) has shown how to compute the PN expansion of the radial integral ( 3.12|) 
to any order in the 1/R expansions 



). At the present 2PN order, Ref. [[LlJ gave a general 
formula (their equation (3.9)) which can be straightforwardly applied to our case. 

As we said above, the function describing the "energy levels", £q = J 7 (No, Jo), is a 
coordinate-invariant construct. As a check on our calculations, we have computed it (or 
rather, we have computed the radial action I r (£q, Jo)) in the two preferred coordinate gauges 
for a spherically symmetric metric: the "Schwarzschild gauge" and the "Isotropic" one. If 
(ai, bi) denote the expansion coefficients ( |3.4| ) in the Schwarzschild gauge (Cs(R) = 1), we 
find (at the 2PN accuracy) 

3/2 



Ir{£o, Jo 



-2£ NR 



<?NR 

A + B 



itlq cr 



C 



pNR 
c 

mo c 2 



-Jo 



+ 



m 



c 2 J 



cNR 



m cr 



+ 



ml 



c*J 3 



where £q R = £q 



mo c 2 , and where 



A 



D 



"2 Ql 



«2 

~2~ 



B 
ai&i 



bi - - ai 

o 



E 



a 1 — a 2 — 







_ b i 


19 


4 


~64 ai ' 


a\bi 


bi h 


2 


8 2 



(3.13) 



64 



[24 a\ - 48 a{ a 2 + 8 a 2 . + 16 a x a 3 - 8 af b% + 8 a 1 a 2 b x - a{b{ + A a{ b 



(3.14) 



Denoting by (a$, fej) the expansion coefficients (|3.4[ ) in the Isotropic gauge (Cj(R) = 
Bi(R)), we find, by calculating directly in the isotropic gauge, that the coefficients 
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A, B, . . . , F entering Eq. ( p.!3| ) have the following (slightly simpler) expressions in terms 
of a,- and : 



A 



1~ 
2 ' 



r h 19 _ 



_ a? a 2 ai&i 



B = bi - - at , 

o 

= a\ — a 2 — ai&i + 6 2 



1 



F = - [3 a x — 6 a 1 a 2 + a 2 + 2 ai 03 — 4 61 + 4 ax a 2 b\ + a 1 b l + 2 a x 6 



(3.15) 



The numerical values of the coefficients A, B, . . . , F are checked to be coordinate-invariant 
by using the following relation between the (a*, bi) and the (ctj, bi) (which is easily derived 
either by integrating dRj/Rj = a/ B s (R s ) dR s /Rs or by using the algebraic link R s = 
Rj y/Bj(RPj) 

a 1 =a 1 , h = b 1 , 



a 2 = a 2 - - ax h , 



b 2 = -b 2 --b 1 



7 u2 1 . 

a 3 = a 3 - a 2 bi + — a 1 b 1 - - ax b 2 . 

lb 4 



(3.16) 



Finally, solving iteratively Eq. ( |3.13| ) with respect to £(f R , we find the analog of Eq. ( |2.15| ), 
i.e. the explicit formula giving the effective "energy levels". It is convenient to write it in 
terms of M = 1^ + J , of the coupling constant 



«o = G M m 



(3.17) 



where M Q is defined by Eq. ( |3.5| ), and of the (GM )-rescaled, dimensionless expansion coef- 
ficients cii and 6j, of the Schwarzschild gauge: 



ai = ai/{GM )\ 



bi/{GM ) 



(3.18) 



with ax = —2. 
We find 



€o(Mo,Jo) = m c - 



2 1 mo a 



2 N* 



1 + 



a ( ^3,1 _|_ C*4,0 



c 2 \MoJ Ml 

Ga2 C§x , Cqq 



(3.19) 



where the coefficients C Pi9 (which parametrize the contributions cx — \ (ao/c) p+q Af ~ p J Q ~ 9 to 
£o/m c 2 ) are given by 
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C*3,l 
C3,3 
G5A 



2D, 
2F, 

-(ABD + E) 



£4,0 — —B , 
C40 = 3D 2 , 



a 



do 



-(SB- 



2C 



(3.20) 



Here, the dimensionless quantities B,C,D,E,F are the GM -rescaled versions of the co- 
efficients of Eq. (|3.13| ), given by replacing the a^s by <Zj in Eqs. fl3.14|) . For instance, 
B = bx - 7/8 ai = 61 + 7/4, etc.. 



IV. RELATING THE "REAL" AND THE "EFFECTIVE" ENERGY LEVELS, 
AND DETERMINING THE EFFECTIVE METRIC 

We still have to define the precise rules by which we wish to relate the real two-body 
problem to the effective one-body one. If we think in quantum terms, there is a natural 
correspondence between Af and Afo, and J and j7~o, which are quantized in units of h. It is 
therefore very natural to require the identification 

Af = Afo, J = Jo, (4.1) 

between the real action variables and the effective ones, and we will do so in the following. 
What is a priori less clear is the relation between the real masses and energies, m 1; m 2 , 
^reai = ( m i + m 2) c 2 + an d the effective ones, mo, M , So = m c 2 + £q R . The usual 
non-relativistic definition of an effective dynamics associated to the relative motion of a 
(Galileo-invariant) two-body system introduces an effective particle whose position q is 
the relative position, q = q 1 — q 2 , whose inertial mass m,Q R is the "reduced" mass fi = 
itii m 2 /(mi+m 2 ), and whose potential energy is the potential energy of the system, V e s(q ) = 
V^eai(<Zi — <7 2 )- m * ne present case of a gravitationally interacting two-body system, with 
KeEd = —G mim 2 /\q 1 — q 2 \, this would determine 

m^ R = /i, and M NR = m x + m 2 = M , (4.2) 

such that a rea i = Gm\ = ao = G MQ R m^ R . The non-relativistic identifications (|4.2j) are, 
however, paradoxical within a relativistic framework, even if they are modified by "relativistic 
corrections", so that, say, m = /x + 0(c~ 2 ), M = M + 0(c~ 2 ), because the reference level 
(and accumulation point for Af, J — > 00) of the real relativistic levels ( 2.15Q will be the 
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total rest-mass-energy Mc 2 , and will therefore be completely different from the reference 
level m Q c 2 ~ /ic 2 of the effective relativistic energy levels. This difference in the relativistic 
reference energy level shows that, while it is very natural to require the straightforward 
identifications ( |4.1|) of the action variables, the mapping between £ real and S must be more 
subtle. 

One might a priori think that the most natural relativistic generalization of the usual 
non-relativistic rules for defining an effective one-body problem consists in requiring that 

£o(A/"o, Jo) = E leal (Af, J) -co, (4.3) 

with a properly chosen constant Cq = M c 2 — mo c 2 taking care of the shift in reference 
level. The rule ( |4.3j ) is equivalent to requiring the identification of the "non-relativistic" 
Hamiltonians (with subtraction of the rest-mass contribution) 

Hr(q', P ') = HZ(q,p), (4-4) 

where the canonical coordinates in each problem must be mapped (because of the identifi- 
cation ( [4.1| )) by a canonical transformation, 

^ d< ? = d ^ + d ^ > ( 4 - 5 ) 

i i 

with some "generating function" g(q,q'). 

We have explored the naive identification (fOf), or (f4.4|) , and found that it was unsatis- 
factory. Indeed, one finds that it is impossible to require simultaneously that: (i) the energy 
levels coincide modulo an overall shift ( |4.3j ), (ii) the effective mass mo coincides with the 
usual reduced mass fi = mi m 2 /(mi + m 2 ), and (iii) the effective metric (|3~3| ) depends only 
on mi and m<i- [This impossibility comes from the fact that the requirement ( f4.4|) is a very 
strong constraint which imposes more equations than unknowns.] If one insists on imposing 
the naive identification ( |4.3| ) there is a price to pay: one must drop at least one of the re- 
quirements (ii) or (iii). Various possibilities are discussed in the Appendices of this paper. 
One possibility is to drop the requirement that mo = [i. As discussed in App. |A], we find 
that there is a unique choice of masses in the effective problem, namely 

m Q = n f-2 ; GMq = GM £ 3 , (4.6) 

with 
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s 5 L 



2V100 + 30z/ + 4z/ 2 - 15 



v 



(4.7) 



which is compatible with the requirements (i) and (iii) above. However, we feel that it is 
quite unnatural to introduce an effective mass mo which differs from /1 even in the non- 
relativistic limit c — > +00. We feel also that this possibility is so constrained that it is only 
available at the 2PN level and will not be generalizable to higher post-Newtonian orders. 
A second (formal) possibility is to introduce some energy dependence, either in mo, say 

pNR /pNR\ 2 



or in the effective metric Namely, the various coefficients ax, bx, (12, b%, a Zi ■ ■ ■ in Eq. 

can be expanded as 

pNR / pNR \ 2 

ai (S ) = a 0) + a?> + a 4) ( ^4 + ■ ■ • , (4.9) 

etc. These possibilities are discussed, for completeness, in App. 0. 

Though the trick of introducing an energy dependence in (both) itlq and the effective 
potential has been advocated, and used, in the quasi-potential approach of Todorov |§, @, 
we feel that it is unsatisfactory. Conceptually, it obscures very much the nature of the map- 
ping between the two problems, and, technically, it renders very difficult the generalization 
(we are interested in) to the case where radiation damping is taken into account (and where 
the energy is no longer conserved). We find much more satisfactory to drop the naive re- 
quirement QjO ), and to replace it by the more general requirement that there exist a certain 
one-to-one mapping between the real energy levels and the effective ones, say 

£o(A/"o, Jo) = f [^cal(AT, J)} . (4.10) 

In explicit , expanded form, the requirement ( 4.10|) yields a deformed version of Eq. (|4.3|) : 



<?NR pNR / <?NR /pNR\ 2 

^0 _^eal ( 1 + ai ^ + Q Jbs^ +■••] . (4.11) 



m c 2 \x c 2 \ [ic 2 \fj,c 2 J 



Here, we assume that the standard identification ( |4.3D holds (together with tuq = fM+0(c~ 2 )) 
in the non-relativistic limit c — » 00. 

We are going to show that the a priori arbitrary function /, i.e. the parameters cti, 0:2, ■ ■ ■ 
can be uniquely selected (at the 2PN level) by imposing the following physically natural 
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requirements: (a) the mass of the effective test particle coincides with the usual reduced 
mass, 

ma = (i, (4.12) 

and, (b) the linearized ( "one-graviton-exchange" ) effective metric coincides with the lin- 
earized Schwarzschild metric with mass M = mi + ma, i.e. 

ai = -2GM , bx = 2GM . (4.13) 



Note that the requirement (|4.12 ) is actually imposed by dimensional analysis as soon as one 
requires m = // + 0(c~ 2 ). Indeed, as we bar any dependence on the energy, it is impossible 
to write any correction terms 0(c~ 2 ) in the link between mo and /i. The requirement ( 4.13Q is 
very natural when one thinks that the role of the effective metric is to reproduce, at all orders 
in the coupling constant G, the interaction generated by exchanging gravitons between two 
masses mi and m^. The "one-graviton-exchange" interaction (linear in G mxm-i) depends 
only on the (Lorentz-invariant) relative velocity and corresponds to a linearized Schwarzschild 
effective metric in the test-mass limit v — > 0. As the coefficient —\a,i is fixed (by dimensional 
analysis, as above) to its non-relativistic value — | ai mo = GMqitiq = Gmim2, it is very 
natural not to deform the coefficient b\ by ^-dependent corrections. 

Let us now prove the consistency of the requirements ( f4.12| ), ( 4.13 ) and determine the 
energy mapping /. We can start from the result ( |3.13j ), in which one replaces Sq R by the 
expansion ( |4.11| ). This leads again to an expression of the form (|3.13|) , with S^ R replaced 

•NR 
real" 



by S^ e R . One can simplify this expression by working with scaled variables 



rO rreal <?NR pNR 

fO _ 1 R frcal _ 1 Ft _ • TP = ~2 P — real 

In — , In — — t r , £/q — ) -C'real — j 

a a m /i 

k = ~, 3 = ~- (4-14) 

«0 Ci 

Here ao = GM Q m and a = GM fi = G nil m 2 as above. We use also the scaled metric 
coefficients and bi of Eq. ( |3.18| ). Let us note, in passing, that, very generally, the di- 
mensionless quantity S /c 2 = £ / (m c 2 ) = 1 + c -2 E is expressible entirely in terms of the 
dimensionless scaled action variables I^/c— /°/(ao c) and of the dimensionless scaled metric 
coefficients a,, 6j. [This scaling behavior can be proved very easily by scaling from the start 
the effective action S = — f rriQcds^ = —a Q c J ds^f with ds^ = (GM )~ 2 ds%, and by using 
scaled coordinates: R = R/GM , t = t/GM .] 
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Let us now make use of the assumptions m = fi and GM = —\a\ = GM (so that 
«o = GMquiq = GM \i = a). But, let us not yet assume the second equation (|4.13|) , i.e. let 
us assume a x = —2, but let us not yet assume any value for b\ = bi/GM = b\/GM. Within 
these assumptions, the scaled version of the result (|3.13|) , with Sq R replaced by ( f4.11|) , reads 

2" 



1 



ll{E {E Ieal ),j ) = -j + 

■ 4t- D - E 



c 2 j 



\/—2E Tea _i 
E re3j \ 



+ 



E- 



real 



c 4 f 



(4.15) 



where 



A = <2i = 1 , 

2 

d , = l9 h a, 
32 4 2 



B' 



7 
4 

7' 
4 



+ 



a 1 — 



«2 



(4.16) 



and where D, E and F are obtained from the expressions ( |3.14| ) by the replacements a« — > a^, 
6i -> &< (with ax = -2). Finally, identifying [7°(£ , ^o)]^^ with I^eai, Jreal), or 
equivalently I^(E (E res i), j ) with I R (E real , jo), yields five equations to be satisfied, namely 
the equations stating that B', C, D, E and F coincide with the corresponding coefficients 
in Eq. ( |2.14| ). The explicit form of these equations is 



7 



4 
19 



a i 

T 



h «i 

32 4 2 

a 2 b\ 
2 - — + — = 
2 2 



15 

T 

b x 



a 2 35 15 3 o 

— = 1 v H f 

2 32 16 32 



4 - a 2 + h 



= 3 



+ 



15 

T 



3z/ . 



6 — 3 a 2 H 

8 2 



bi - -a 2 h 



6? b_2 

16 4 



35 

T 



V . 



(4.17) 
(4.18) 

(4.19) 
(4.20) 
(4.21) 



Note that the subsystem made of the two equations ( |4.17j) , ( |4.18| ) (corresponding to B' 
and C) contains the three unknowns b±, a±, a 2} while the three equations ( |4.19| )- (|4.21 



(corresponding to D, E and F) contains the unknowns b\, b 2 , a 2 , a^. In this section we shall 
consider only the first ( il BC" ) subsystem, leaving the U DEF" system to the next section. 

It is easily seen that the BC subsystem would admit no solution in b\ if we were to 
impose ot\ = a 2 = 0. This proves the assertion made above that one needs a non trivial 
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energy mapping £ = /(£ rea i). On the other hand, if we introduce the two free parameters ol\, 
a 2 the BC subsystem becomes an indeterminate system of two equations for three unknowns. 
As argued above, it is physically very natural to impose that the linearized effective metric 
coincides with the linearized Schwarzschild metric, i.e. that 



6i = 2. (4.22) 
Then the BC system ( |4.17| ), ( |4.18|) admits the unique solution: 



2 

This solution corresponds to the link 



ai = -, a 2 = 0. (4.23) 



<?NR <?NR / „ <?NR\ 
°0 °real ' - u °™ 1 > 



which is equivalent to 



m c 2 /ic 2 \ 2 [i c 2 J ' 

£0 = gaj -mjc 4 - m 2 c 4 
mo c 2 2 mi m 2 c 4 



1 + - , (4.24) 



(4.25) 



Remarkably, the map ( |4.25| ) between the real total relativistic energy S Tea \ = M c 2 + 



NR 
real' 

and the effective relativistic energy £ = m Q c 2 + £q R coincides with the one introduced 
by Brezin, Itzykson and Zinn- Justin which maps very simply the one-body relativistic 
Balmer formula onto the two-body one (in quantum electrodynamics). The same map was 
also recently used by Damour, Iyer and Sathyaprakash ||21|| . There it was emphasized that 
the function ip(s) of the Mandelstam invariant s = S 2 eal appearing on the R.H.S. of Eq. ( |4.25| ) 
is the most natural symmetric function of the asymptotic[] 4-momenta p±, p% of a two-particle 
system which reduces, in the test-mass limit m 2 mi, to the energy of m 2 in the rest-frame 
of mi. Indeed, (setting here c = 1 for simplicity) 

/ j _ s - m i - m l _ + Vif - ml - ml _ pi-p 2 
2m±m 2 2m\m 2 m\m 2 

Finally, we have two a priori independent motivations for using the function (p(s), i.e. the link 
( |4.25| ), to map the real two-body energy onto the effective one-body one: (i) the simplicity, 



1 We consider here scattering states. By analytic continuation in s, the function (f(s) is naturally 
expected to play a special role in the energetics of two-body bound states. 
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and the symmetry, of the expression (|4.26|) which generalizes the test-mass conserved energy 
So/mo = —K^pQ/m (where is the Killing vector defined by the time direction of the 
background field) (see |21|]), and (ii) the fact that it corresponds to a linearized effective 
metric coinciding with the linearized Schwarzschild metric. Actually, these two facts are 
not really independent, because (as discussed in []IJ and 0) they correspond heuristically to 
saying that the "effective interaction" is the interaction felt by any of the two particles in 
the rest frame of the other particle. 

Summarizing: The rules we shall assume for relating the real two-body problem to the 
effective one-body one are Eqs. fl4.1|) (or equivalently the condition ( |4.5|) that the phase-space 
coordinates are canonically mapped), and Eq. (|4.25| ). 



V. THE EFFECTIVE ONE-BODY METRIC AND THE DYNAMICS IT DEFINES 

Having specified the rules linking the real two-body problem to the effective one-body 
one, we can now proceed to the determination of the effective metric (at the 2PN level). We 
shall work in Schwarzschild coordinates: 

ds 2 eS = -A(R) c 2 dt 2 + B(R) dR 2 + R 2 (d6 2 + sin 2 9 dtp 2 ) , (5.1) 

with A(R) and B(R) constructed as expansions of the form (|3.4j ). It will be useful to rewrite 



also the effective metric in the form 

ds 2 eS = -A(R) c 2 dt 2 + dR 2 + R 2 (d6 2 + sin 2 9 dip 2 ) , (5.2) 

A(R) 

in which we factorize, a la Schwarzschild, in front of the dR 2 term, and consider that, 
besides A(R), the second function constructed as an expansion in 1/R is 

D(R) = A(R) B(R) = 1 + + -^L + . . . , (5.3) 

where 

di — di + 61 , d 2 = a 2 + a x hi + b 2 . (5.4) 

To determine the effective metric, i.e. the coefficients aj and bi, or equivalently and 
di = di/(GM) 1 , we insert the known values of 61, a± and a 2 (namely b\ — 2, a,\ — is/2, 
a 2 = 0) into the remaining equations ([4.19|) - ([4.21|) ( U DEF system"). This yields three 
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equations for the three unknowns 0,2, hi and 03. The unique solution of this DEF system 
reads 

a 2 = , a 3 = 2u , b 2 = 4 - 6z/ . (5.5) 

In other words, our natural assumptions ( |4.12|) , ( [4.1 3|) have led us uniquely to the simple 



energy map (|4.25|) and to an effective one-body metric given by 

3 



2GM ( GM 
2GM ,. _ , / GM\ 2 



B (*) = l + _ + (4-6„)(^j +..., (5.7) 
D(R) = l-sJf^J+--. (5.8) 



The simplicity of the final results ( |5.6| )-( |5T8| ) is striking. The effective metric ( |5.2j ) is a 
simple deformation of the Schwarzschild metric (A S (R) = 1 — 2GM/c 2 R, D S (R) = 1) with 
deformation parameter v. Note also that there are no //-dependent corrections to A(R) at 
the 1PN level, i.e. no u(GM/c 2 R) 2 contribution to A(R). The first //-dependent corrections 
enter at the 2PN level. Remembering that the (2PN) effective metric fully encodes the 
information contained in the complicated 2PN expressions ( 2.1 4j) or ( 2.15| ), it is remarkable 



that the metric coefficients ( ^.6|) -( pl8|) be so simple. The previous approach of Ref. @] led 
to much more complicated expressions at the 1PN level (to which it was limited). 

In this paper, we propose to trust the physical consequences of the effective metric ( |5.2| ) . 
with A(R) given by Eq. ( |5.6| ) and D(R) given by Eq. ( |5.8| ), even in the region where R 
is of order of a few times GM/c 2 . Note that even in the extreme case where u — 1/4 
and R ~ 2 GM/c 2 the //-dependent additional terms entering the effective metric remain 
relatively small: Indeed, in this case, 5 l/ A(R) = 2v[GM/c 2 R) 3 = 1/16 and —5 l/ D(R) = 
6z/ (GM/c 2 R) 2 = 3/8. We expect, therefore, that it should be a fortiori possible to trust the 
predictions of the effective metric ( |5.2D near the innermost stable circular orbit, i.e. around 
R ~ 6 GM/c 2 (where S V A(R) - 2x 1(T 3 and -8 V D(R) ~ 4 x 1(T 2 ). Note that this nice 
feature of having only a small deformation of Schwarzschild, even when v = 1/4, is not 
shared by the "hybrid" approach of Kidder, Will and Wiseman |22[. Indeed, as emphasized 
in Ref. |ZT| , the //-deformations considered in the hybrid approach are, for some coefficients, 



larger than unity when v = 1/4. This is related to the fact pointed out by Schafer and Wex 
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|23|j , | 24j that, by applying the hybrid approach of |22[ to the Hamiltonian, instead of the 
equations of motion, one gets significantly different predictions. 

Let us note also that, if we decide to write the effective metric in the form ( |5.2| ), the 
existence of a simple zero in the function A(R), say A(R H ) = 0, implies (if D(R H ) ^ 0, 
and D(R) > for R > R H ) that the hypersurface R = Rh is (like in the undeformed 
Schwarzschild case) a regular (Killing) horizon. As usual, one can define Kruskal-like coor- 
dinates to see explicitly the regular nature of the horizon R = Rh (made of two intersecting 
null hypersurfaces) . In our case, one checks easily that the function v42pn(-R) defined by 
the first three terms on the R.H.S. of (|5.6| ) admits a simple zeroQ at some Rh{v), when 
< v < |. The position Rr{v) of this "effective horizon" smoothly, and monotonically, 
evolves with the deformation parameter v between Rh{0) = 2GM/c 2 and 

R H (1/4) ~ 0.9277 (2GM/c 2 ) . (5.9) 

This relatively small change of the horizon toward a smaller value, i.e. a smaller horizon area 
(to quote an invariant measure of the location of the horizon), suggests that the dynamics 
of trajectories in the effective metric will also be only a small deformation of the standard 
Schwarzschild case. 

One of the main aims of the present work is indeed to study the dynamics (and the 
energetics) in the effective metric ( |5.2|) . In particular, as gravitational radiation damping is 
known to circularize binary orbits, we are especially interested in studying the stable circular 
orbits in the effective metric. A convenient tool for doing this is to introduce an effective 



potential [pq], p9|. Note that the Hamilton- Jacobi equation (|3.9|) yields 



£ ° ^-Wj (R) + ^(^) >W Jo (R), (5.10) 



m^c 2 J B(R) \m c y 

where Pr = dS c s/dR is the effective radial momentum, and where the "effective radial 
potential" Wj (R) is defined as 

(Jo/ mo c) 2 " 



W Jo (R)=A(R) 



1 + 



C(R) R 2 



(5.11) 



2 We consider here only the zero of ^2PN (R) which is continuously connected to the usual horizon 
R S H = 2GM/c 2 when v -> 0. 
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We read also from Eq. (|5.10|) the relativistic effective Hamiltonian 



H*(R,P R ,P V 



m c 




+ 



p 2 



2 c 2 B(R) m 2 c 2 C(R)R 2 



= m c 2 Jw P (R) + 



A(R) ( P R 



B{R) \moC 



(5.12) 



The coordinate angular frequency along circular orbits is obtained by differentiating the 
Hamiltonian, that is 



\dt), 



( dH*{R,P R ,P v 
V dP v 



(5.13) 



Pr=0 



which gives explicitly (using P^ = J ) 



UJ 



Jo 



m C(R) R 2 



1 + ss 



■j 2 



(5.14) 



C{R) R 2 



Eqs. (|5.11| ) and ( 5.14 ) are valid in an arbitrary radial coordinate gauge, but we shall use 
them in the Schwarzschild gauge where the metric coefficient C(R) = 1. Note that W(R) 
and ujo then depend only on the metric coefficient A(R). In dimensionless scaled variables 
R = c 2 R/(GM), j = cJ /(GM fi),Q Q = GMlu /c 3 (in our case M = M and m = //), the 
effective potential and the orbital frequency (along circular orbits) are quite simple: 



W Jo (R) = A(R) 



3o_ 
R 2 



R 2 . 

A(R) 



i + S- 

^ R 2 



(5.15) 



If we define the 2PN-accurate A(R) by the straightforward truncation of Eq. (|5.6|) , namely 

A 2PN (R) = l-l + % 1 (5.16) 



Wj is a fifth-order polynomial in u = l/R = GM/(c 2 R). As the analytical study of 
the extrema of Wj is rather complicated, we have used a numerical approach. When v 
varies between and 1/4, Wj evolves into a smoothly deformed version of the standard 
Schwarzschild effective potential. To illustrate this fact, we plot in Fig. [l| Wj (R) for v = \ 
and for various values of the dimensionless angular momentum j Q . Note that the latter 
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v = 1/4 




30.0 



FIG. 1. The effective radial potential Wj(R) (at the 2PN level and for v = 1/4) versus the 
dimensionless radial variable c 2 R/(GM) for three different values of the dimensionless angular 
momentum j = cj xea \j (GM [i) . Note that the effective radial potential tends to one for R — > oo. 
The stable circular orbits are located at the minima of the effective potential and are indicated by 
heavy black circles. The innermost stable circular orbit corresponds to the critical value j*. In the 
case of the j\ curve the orbit of a particle with energy Eq = £q is an elliptical rosette. 

quantity coincides (in view of our rules) with the corresponding real two-body dimensionless 
angular momentum j: 

G M m G M fi 

[Note that our definition of the j's differs by a factor c from the one used in the previous 
section.] 

As usual, because of the inequality ( |5.10| ), when j and So = 80/(1110 c 2 ) are fixed, the 
trajectory of a particle following a geodesic in the effective metric fl5.2|) can be qualitatively 
read on Fig. |l|. For instance, in the case illustrated for the j\ curve (E* = So line), the 
orbit will be an elliptical rosette, with the radial variable oscillating between a minimum 
and a maximum (solid line in Fig. [I]). The stable circular orbits are located at the minima 
of the effective potential (the maxima being unstable circular orbits). The innermost stable 
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circular orbit (ISCO) corresponds to the critical value j t of the angular momentum where 
the maximum and the minimum of the effective potential fuse together to form an horizontal 
inflection point: 







d 2 W. 



am 



(5.1* 



Let us, for comparison with our deformed case, recall the standard results for circular orbits 
in a Schwarzschild spacetime [28], [ 29| . With the notation u = GM /c 2 R (for a Schwarzschild 



metric of mass M ), the location, orbital frequency^, and energy of circular orbits are given, 
when j varies, by 



u 



1 - Wl 



12 



J 



(5.19) 



~ GM 3/2 
uj s = — =- u = u 1 , 

c 6 



(5.20) 



En 



mo c z 



j(l-2u)u 



1/2 



The ISCO corresponds to the critical values 

c 1 



J * 



1 



(5.21) 



(5.22) 



6' ~* 6^6' ~* V9 
In the deformed Schwarzschild case defined by Eq. ( |5.16|) , the ISCO for the extreme case 
v = ~ is numerically found to correspond to the values 



1 2PN 



,2PN 



U * = 



2PN 

0* 



?2PN 
c 0* 



GM^Jisco 
f GM\ 

\c R) ISCO 
V c 3 



ISCO 



= 3.404 = 0.983 j; 
0.1749 = 1.049 wf , 
= 0.07230 = 1.063 Of, 

0.94040 = 0.99744 £f. 



A 4 c / isco 



(5.23) 
(5.24) 
(5.25) 
(5.26) 



3 Here, as well as in Eqs. ( 5.2 5| ) and Q5.31 ) below, uj denotes the angular frequency dtp/dt on a 
circular orbit (in the equatorial plane). 
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Note that the Schwarzschild-coordinate radius of the effective ISCO is (when v = 1/4) 
R lsco = 5.718 GM/c 2 , i.e. lower than the standard Schwarzschild value 6GM/c 2 corre- 
sponding to the total mass M = ni\ + 777.2. This is consistent with the fact that the effective 
horizon was drawn in below 2GM/c 2 when v was turned on. Note, however, that the three 
quantities Mq PN , ^0™ an d £0™ entering equations (|5.24|) - (|5.26|) are mathematical quanti- 
ties defined in the effective problem, and not physical quantities defined in the real problem 
(hence the subscript added as a warning). [By contrast, j 2PN , Eq. ( |5.23| ) is directly related 
to the real, two-body angular momentum.] For physical (and astrophysical) purposes, we 
need to transform the information contained in Eqs. ( |5.24|) -( |5.26|) into numbers concerning 
physical quantities defined in the real, two-body problem. For the energy, this is achieved (by 
definition) by using Eq. fl4.25|) to compute the real, two-body total energy £ Tea \. Explicitly, 
the solution of Eq. ( |4.25|) is (see also 0) 



£eal = MC 2 \l 1 + 2 »( S0 m ™f )- 

We need also to transform the effective orbital frequency ujq. This is easily done as follows. 
We know that the Hamiltonians of the real and effective problems are related by a mapping 

iWC al ) = M#o(4°)) , (5-28) 

where a = R,9,ip (for the 3-dimensional problem), and where the function h (the inverse of 
the function / of Eq. ( f4.10| )) is, in our case, explicitly defined by Eq. ( [5.271 ). On the other 
hand, we know that the action variables are identically mapped onto each other: J° = 7^ eal 
(canonical transformation). The frequency of the motion of any separated degree of freedom 
is given by the general formulas cu® = d H (I°)/d I®, uo T ™ x = d H rcal (r eBl ) / d F a e£L \ where the 
Hamiltonians are considered as functions of the canonically conjugate action- angle variables 
(I a , 9 a) (remembering that for such integrable systems, the Hamiltonian does not depend on 
the 9 , s). Therefore the frequencies of the real problem are all obtained from the frequencies 
of the effective one by a common, energy- dependent factor 

u?* dtp dH rcal _ dh(H ) 

to° a dt^ 1 dH dH " 1 ' 

In our case this "blue shift" factor reads 



4 For bound states, k/ eal > co°. 
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FIG. 2. Variation with v (at the 2PN level) of the ISCO values of the real non-relativistic 
energy E Tea \ = £^ = (£ rea i — Mc 2 )//xc 2 (on the left) and of the real angular momen- 
tum j = cJ vea \/GMn. (on the right), divided by the corresponding Schwarzschild values 
\Es\ = \£^ R \ = 1 - ~ 0.05719 and j s = Vl2, respectively. 



real 



dtp, 



(5.30) 



"2 dt^ y /l + 2is(£ -m c 2 )/m c 2 ' 

As indicated in Eqs. ( |5.29| ) and ( |5.30| ) the same energy-dependent "blue shift" factor maps 
the effective and the real times (along corresponding orbits). Note that we have here a 
simple generalization of the spatial canonical transformation (dp A dq = dp A dq ) to the 
time domain (dH A dt = dHo A dto). 

Applying the transformations ( |5.27| ) and ( |5.29| ), we obtain the physical quantities^ pre- 
dicted by our effective 2PN metric, still in the extreme case v = 1/4, 



5 In Eq. (|5T31p ^reai — dE xc& \j dj xe& \ is again the angular frequency on a circular orbit. It should 
not be confused with the radial (periastron to periastron) frequency lor for non-circular, rosette 
orbits. 
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FIG. 3. ISCO values (at the 2PN level) of the quantity z = (GMw^/c 3 )" 2 / 3 , divided by the 
Schwarzschild value zs = 6, versus v. 



to. 



-2PN 
real* 



(GM 



^rcal 



1.079 ZD? = 0.07340, 



ISCO 



1.050 (5; 



1' 



-0.06005 



(5.31) 
(5.32) 



^ C / ISCO 

We represent in Figs. § and |3] the variation with v of the ISCO values of the real non- 
relativistic energy, E re& \ = = (^ rca i — M c 2 )/jj,c 2 } the real angular momentum, j = 
cJ^reai/GMfi, and of the quantity 

-2/3 

(5.33) 



GM 



V c 3 



■ ^real 



which is an invariant measure of the radial position of the orbit, and which coincides with 
the scaled Schwarzschild radius R = c 2 R/(GM) in the test-mass limit v — > 0. One checks 
that our ISCO values respect the "black hole limit" J^aX < GS 2 eal /c 5 , so that the system 
does not need to radiate a lot of gravitational waves in the final coalescence before being 
able to settle down black hole. 

Let us now briefly compare our predictions with previous ones in the literature. The 
first attempt to address the question of the ISCO for binary systems of comparable masses 
was made by Clark and Eardley |30]]. They worked only at the 1PN level, and predicted 
that the ISCO should be significantly more tightly bound than in the Schwarzschild case 
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(with M = M = mi + m 2 ): Sq^/ ficr ~ —0.1 when u — 1/4, compared to £schwarz/ m o c 



a/8/9 — 1 ~ —0.0572. Blackburn and Detweiler [31] used an initial value formalism (which 
is only a rough approximation, even in the test-mass limit) to predict an extremely tight 
ISCO when v = 1/4 : S^/fic 2 ~ —0.7. Kidder, Will and Wiseman Q were the first to try 
to use the full 2PN information contained in the Damour-Deruelle equations of motion (|1 . 1| ) 
to estimate analytically the change of the ISCO brought by turning on a finite mass ratio 
v. They introduced an "hybrid" approach in which one re-sums exactly the "Schwarzschild" 
(^-independent) terms in the equations of motion, and treats the //-dependent terms as 
additional corrections. In contrast with our present 2PN-effective approach (and also with 
the less reliable previous studies [0O[| , |3lj), they predict^ that, when v increases, the ISCO 



becomes markedly less tightly bound: e.g. £^ w //ic 2 ~ —0.0377 when v = 1/4. If their 
trend were real, this would imply that, except for the very stiff equations of state of nuclear 
matter (leading to large neutron star radii), the final plunge triggered when the ISCO is 
reached by an inspiraling (1.4M + 1.4M ) neutron star binary would probably take place 
before tidal disruption. However, both the robustness and the consistency of the hybrid 



approach of |22[] have been questioned. Wex and Schafer |23| showed that the predictions 
of the hybrid approach were not "robust" in that they could be significantly modified by 
applying this approach to the Hamiltonian, rather than to the equations of motion. Schafer 



and Wex [24] further showed that the predictions of the hybrid approach were not robust 



under a change of coordinate system. Moreover, Ref. |21[] has questioned the consistency of 
the hybrid approach by pointing out that the formal 'V-corrections" represent, in several 
cases, a very large (larger than 100%) modification of the corresponding //-independent terms. 
This unreliability of the hybrid approach casts a doubt on the ISCO estimates of Ref. p5 



which are based on hybrid orbital terms, and which use only 1PN accuracy in most terms. 

Damour, Iyer and Sathyaprakash [^1[] have introduced (at the 2PN level) another ana- 
lytical approach to the determination of the ISCO, based on the Pade approximants of some 
invariant energy function (closely related with the energy transformation ( |4.25| )). Their 



trend is consistent with the one found in the present paper, namely a more tightly bound 



We use here the values read on the figures 3 and 4 of Ref. ||: for £ NR and (mf) = 0.00963, 
which refer to a static ISCO without radiation damping. 
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FIG. 4. ISCO values (for v = 1/4 ) of the real non-relativistic energy 

E = £j£*> divided by 

the corresponding Schwarzschild value Es = £3^, versus z/zs- On the left we have compared our 
predictions at the 1PN level (■) and 2PN level (♦) with the results obtained in [21] (►) and [22] (m). 
The (*) indicates the Schwarzschild predictions. The right panel is a magnification of the part of 
the left one in which we analyze the robustness of our method by exhibiting the points (•) obtained 
by introducing in the effective metric reasonable 3PN and 4PN contributions: (04, a' 5 ) = (±4, —4), 



(±4,0) and (±4, +4) in the notation of Eq. (5.34). 



ISCO: for v = 1/4, the Pade approximant approach predicts S^/fic 2 ~ —0.0653. 

Numerical methods have recently been used to try to locate the ISCO for binary neutron 
stars |p6|| , f27| . However, we do not think that the truncation of Einstein's field equations (to a 
conformally flat spatial metric) used in these works is a good approximation for close orbits. 
Indeed, at the 2PN approximation, some numerically significant terms in the interaction 
potential come from the transverse-traceless part of the metric [13], 0, [10]. Moreover, 
the (unrealistic) assumption used in these works that the stars are corotating has probably 
also a significant effect on the location of the ISCO by adding both spin-orbit and spin-spin 
interaction terms. 

This large scatter in the predictions for the location of the ISCO for comparable masses 
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Method 






z 




/© (kHz) 


"Schwarzschild" 




-0.01430 


6 


0.06804 


2.199 


Eff. action 1PN 




-0.01440 


5.942 


0.06904 


2.231 


Eff. action 2PN 




-0.01501 


5.704 


0.07340 


2.372 


Eff. action (04,05) 


= (-4,-4) 


-0.01462 


5.891 


0.06994 


2.260 


Eff. action (a 4 ,o 5 ) 


= (-4,0) 


-0.01469 


5.854 


0.07061 


2.267 


Eff. action (04,05) 




-0.01476 


5.815 


0.07131 


2.304 


Eff. action (a 4 , a 5 ) 


= (+4,-4) 


-0.01530 


5.583 


0.07582 


2.450 


Eff. action (a 4 , a 5 ) 


= (+4,0) 


-0.01540 


5.531 


0.07688 


2.484 


Eff. action (a 4 , a 5 ) 


= (+4, +4) 


-0.01551 


5.475 


0.07806 


2.522 


D.I.S. [21] 




-0.01633 


5.036 


0.08850 


2.860 


K.W.W. [22] 




-0.00943 


6.49 


0.0605 


1.96 



TABLE I. Summary of the ISCO values used in Fig. 4 (u = 1/4). Note that we give here 
S^^/Mc 2 , that is the ratio between the energy that can be radiated in gravitational waves before the 
final plunge and the total mass-energy initially available. The first row refers to the naive estimate 
defined by a test particle of mass fi in a Schwarzschild spacetime of mass M. We show also in the 
last column the solar-mass-scaled orbital frequency f & defined by / rca i = uj Tca ,\/ (2ir) = f & (M Q /M). 



poses the question of the "robustness" of our new, effective-action approach. The main 
problem can be formulated as follows. Assuming that the effective-action approach (for 
the time-symmetric part of the dynamics) makes sense at higher post-Newtonian levels, the 
"exact" effective function A(R) will read 

The question is then to know how sensitive is the location of the ISCO to the values of the 

(still unknown) coefficients a 4 , a' 5 , One should have some a priori idea of the reasonable 

range of values of a 4 , a' 5 , . . .. A rationale for deciding upon the reasonable values of a' A is the 
following. At the 2PN level, it is formally equivalent to use (with u = GMj c 2 R) A 2 pn = 
l—2u+2uu 3 or the factorized form A' 2PN = (1— 2u)(l+2uu 3 ). However, A' 2PN = A 2 pn~ 4:uu 4 
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which corresponds to a' A = —4. This suggest that —4 < a' A < +4 is a reasonable range. We 
shall also consider —4 < a' b < +4 as a plausible range. Note that both choices correspond to 
having coefficients of u n which vary between —1 and +1 when v = 1/4. The robustness of 
our effective-action predictions against the introduction of a' A and a' 5 is illustrated in Fig. ||| 
The numerical values used in Fig. [| are exhibited in Tab. |. 

Fig. | plots the ratio E/\E S \ where E = S^/fic 2 = (£ rcal - Mc 2 )//2C 2 at the ISCO (for 
v = 1/4) and Eg = a/(8/9) — 1 ~ —0.05719 is the corresponding "Schwarzschild" value, 
versus z/zs where z is defined in Eq. ( (5.33Q , and where zs = 6. This figure compares the 



predictions of Ref . [22 j , of Ref . |2T|] and of our new, effective-action prediction (at the 2PN 
level) . We have also added what would be the prediction of the effective-action approach at 
the 1PN level. Note that, at the 1PN level, the function A(K), Eq. ( |5 . 6| ) , exactly coincides 
with the Schwarzschild one, but that the energy mapping ( |4.24j ) introduces a slight deviation 
from the test-mass limit. Fig. [| exhibits also the points obtained when considering (04, a' 5 ) = 
(±4, —4), (±4, 0) and (±4, +4). We see on this figure that the main prediction of the present 
approach (a prediction already clear from the fact that the 2PN contribution to A(R) is 
fractionally small), namely that the ISCO is only slightly more bound than in the test-mass 
limit, is robust under the addition of higher PN contributions. The sensitivity to a' 4 of 
the binding energy is only at the ~ 3% level (for a' 4 = ±4), while its sensitivity to the 
4PN-coefficient a' 5 is further reduced to the ~ 0.6% level (for a' 5 = ±4). Still, it would 
be important to determine the 3PN coefficient a' A to refine the determination of the ISCO 
quantities. 



VI. EXPLICIT MAPPING BETWEEN THE REAL PROBLEM AND THE 

EFFECTIVE ONE 

The basic idea of the effective one-body approach is to map the complicated and badly- 
convergent PN-expansion of the dynamics of a two-body system onto a simpler auxiliary 
one-body problem. We have shown in the previous sections that by imposing some simple, 
coordinate-invariant requirements, we could uniquely determine that the one-body dynamics 
was defined (at the 2PN level) by geodesic motion in a certain deformed Schwarzschild 
spacetime. The latter dynamics can be solved exactly by means of quadratures (e.g. by using 
the Hamilton- Jacobi method, see Eqs. (|3.7|)-( |3T2"| )). Note that this exact solution defines a 
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particular re-summation of the original 2PN-expanded dynamics. The hope (that we tried 
to subtantiate in Sec. [V]) is that this re-summation captures, with sufficient approximation, 
the crucial non-perturbative aspects of the two-body dynamics, such as the existence of an 
ISCO. 

As all the current work about the equations of motion, and/or the gravitational- wave ra- 
diation, of binary systems is done in some specific coordinate systems (harmonic or ADM), 
we need to complete the (coordinate-invariant) work done in the previous sections by explic- 
itly constructing the transformation which maps the variables entering the effective problem 
onto those of the real one. We have already mentioned that the transformation between 



harmonic and ADM coordinates has been explicitly worked out in Refs. |10| and ||11|| . Here, 
we shall explicitly relate the ADM phase-space variables Q = q l — q 2 and P = dS/dQ 
of the relative motion (as defined in Sec. [TT] above) to the coordinate and momenta of the 
effective problem. More precisely, we shall construct the map 

q H = QV, Pi), Pi = Vitf,Pi), (6.1) 

transforming the reduced ADM relative position and momenta (q l ,Pi), defined in Eq. ( |2.4| ), 
into the corresponding reduced cartesian-like position and momenta {q h ,p'i) canonically de- 



fined by the (Schwarzschild-gauge) effective action ( |3.2| ). In other words, 



O n P' 

<-=§M' * = 7' (6 ' 2) 

with Q' 1 = R sin9 cosp, Q' 2 = R sm9 sincp, Q' 3 = R cos 9, and P[ = dS^/dQ' 1 . Here, 
the "effective" coordinates R, 9, tp are those of Eq. (|5TT| ) (in Schwarzschild gauge) and S c s = 
— J ficds c ff. The corresponding effective Hamiltonian (with respect to the coordinate time t 
of the effective problem) is easily found by solving g^(Q') P'^ P' v + m\ c 2 = in terms of the 
energy £ = — Pq. Transforming the usual polar-coordinate result (equivalent to Eq. ( |5.10| )) 
into cartesian coordinates leads to 



H eS (Q',P')=fic 2 JA(Q') 



(n' ■ P') 2 (n' x P') 2 



fi 2 c 2 B(Q') fi 2 c 2 



(6.3) 



where Q' = 5ij Q h Q'i = R, where n H = Q n /Q' is the unit vector in the radial direction, 
and where the scalar and vector products are performed as in Euclidean space. When scaling 
the effective coordinates as in ( |6.2| ), we need to scale correspondingly the time variable, the 
Hamiltonian and the action of the effective problem: 



32 



Note that the effective Hamiltonian ( |6.3|) contains the rest-mass contribution. The scaled 
version of ( |6.3|) simplifies to 



H eS (q , ,p') = c\ A(q>) 



pi i ■ p'y ( 1 



1+ c 2 + c 2 \B{q>) 1 



(6.5) 



where q' = \/^Jq^7fi = R/GM and n n = q n /q'. As was mentioned above the identification 
of the action variables in the real and effective problems guarantees that the two problems are 
mapped by a canonical transformation, i.e. a transformation such that Eq. fl4.5| ) is satisfied. 
It will be more convenient to replace the generating function g(q, q') of Eq. ( |4.5| ) by the new 
generating function G(q,p') = g(q,q') +v'i1 H such that 

Pl dq l + q' l dp' i = dG(q,p'). (6.6) 

We can further separate G(q,p') into G{&(q,p') = q l p'i, which generates the identity trans- 
formation, and an additional (perturbative) contribution G(q,p'): 

G(q,p') = q l p' l + G(q,p') , G(q,p') = \ G 1PN (q,p') + i G 2PN (Q,p') ■ (6-7) 



& c 4 



Eqs. (KB), (O) yield the link 



H j , dG(q,p') , dG(q,p') 

q=q+ ^r> — a*-- (6 - 8) 

Note that Eqs. (|6.8|) are exact and determine q' and p in function of q and p'. We have, 
however, written them in a form appropriate for determining, by successive iteration, q' 
and p' in function of q and p. If needed (e.g. for applications of the present work to the 
direct numerical calculation of the effective dynamics in the original q,p coordinates), it is 
numerically fast to iterate Eqs. ( |6.8|) to get Eqs. ( |6.1|) . For our present purpose we need an 



explicit analytical approximation of Eqs. (|6.1| ) at the 2PN level. Remembering that G starts 
at order 1/c 2 , one easily finds that 

f_i , dG(g,p) dG(g,p) d 2 G(q,p) | Q f 1 
dpi d qi dpjdpi \c 6 

dG(q,p) dG(q,p) d*G(g,p) (1 
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In the terms linear in G(q,p) one needs to use the full (1PN + 2PN) expression of G(q,p), 
while in the quadratic terms it is enough to use Gipn/c 2 . 

To determine the generating function G(q,p) we need to write the equation stating that, 
under the canonical transformation ( |6.8|) , the effective Hamiltonian H e ff(q',p') is mapped 
into a function of q and p which is linked to the real (relativistic) Hamiltonian H^ al (q,p) 
by our rule (|4.25|) . If we write this link in terms of the reduced effective Hamiltonian (|6.5|) , 
and of the reduced, non-relativistic real Hamiltonian = (-^reai — Mc 2 )/fi (the same as 
H appearing in Eqs. ( |2.5|) , ( |2.6| ) above), it reads 



, + M + » = JL S MM ■ (6.10) 



Actually, we found more convenient to work with the square of Eq. ( |6.1(J| ), so as to get 
rid of the square root in H e g, Eq. (|6.5| ). Hence, writing (half) the square of Eq. ( |6.10| ), 



and Taylor-expanding H e ff[q'(q,p),p'(q,p)} using Eqs. ( |6.7| )-([6J|), we get at order 1/c 4 , the 
following partial differential equation for Gipn (?,£>) 

dH Ncwt dGj PN dH Ncwt 0G 1PN _ v 4 _ p 2 / u\ (n-p) 2 / v\ 1 

dq 1 d Pi d Pi dq i ~ 2 P 1 ] q \ 2J q \ 2) q 2 ' 

(6.11) 

where we have denoted the Newtonian Hamiltonian as H^ ewt = Hq — p 2 /2 — 1/q (see 
Eq. ( |2.6a| )). At order 1/c 6 , a more complex calculation gives the partial differential equation 
for G 2 p N (g,p) 

d-^Newt 9G 2PN dH Ncwt dG 2PN _ v - 3 - - - (n ; p) 2 
—5— ^ ^ ^— — - 77 + I 1 + v ) H o ti2 + Mi - (2 + 3 v) — 

oq l a pi a pi oq l 2 

V 81Z 9G*ipn 97^ 9G*ipn 9Gipn 9 2 Gipn d H^ ewt <9Gipn 9 2 Gipn <9 if Newt 



q 3 dq % dpi dpi dq l dqi dpjdpi dq l dqi dpjdq 1 dp 

1 9 GlPN 9 GlPN 9 2 i?Newt 1 <9 GlPN 9 GlPN <9 2 -f^Newt 



2 <9p,- dqt-dqi 2 <9<f <9g-? dpidpj 



(6.12) 



where if 2 and i?4 are given by Eqs. ( |2.6b|) , ( |2.6c| ), while 

1 



11= - {{n ■ pf + p 2 ) . (6.13) 

q 



The partial differential equations (|6.11|) and ( |6 . 1 2|) have the general form 
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d i^Ncwt 9G n d H^ewt dG n q % dG Ti 



dG r 



Pi 



(6.14) 



dq % dpi dpi dq l q 3 dpi dqi 
where, at each PN order n = 1PN or 2PN, the R.H.S. is a known source term K n (q,p). Note 
that the L.H.S. of Eq. ( |6.14| ) is the Poisson bracket {H^ ewt , G n }, or, equivalently, minus the 
time derivative of G n along the Newtonian motion. It is easily checked that the solution 
of Eq. (|6.14|) is unique modulo the addition of terms generating a constant time shift or a 
spatial rotation. [Indeed, the homogeneous scalar solutions of Eq. (|6.14j ) must correspond 
to the scalar constants of motion of the Keplerian motion: H^ cwt (q,p) and (q x p) 2 .] If 
we require (as we can) that G(q,p) changes sign when q or (separately) p change sign, the 
generating function is uniquely fixed. In particular, at 1PN level, by looking at the structure 
of the source terms, i.e. the R.H.S. of Eq. ( |6.11|) , we can prove in advance that Gipn must 
be of the form 



Gipn(<?,p) = {q-p) 



2 , Pi 
aiP H 

Q . 



(6.15) 



Inserting Eq. ( |6.15| ) in the equation to be satisfied ( |6.11| ) gives a system of four equations 
for the two unknown coefficients ol\ and (3\. Two of these equations give directly the values 
«i and j3i, 



a i 



v 
2 



01 



v 

1 + - 
2 



while the two redundant equations, 

#L = -1 



V 



2«!+/?! = l-- 



(6.16) 



(6.17) 



are identically satisfied by the solution (|6.16| ). 

Using these lPN-results we can go further and evaluate the 2PN-source term K 2 (q,p) in 



Eq. ( grp : 



/<' 2 ( 9 ,p) = ^(l + 3W P 6 +^(-l + 8,)^-^(9 + ,)^^ + |M8 + 3,)^- 
8 8 q 4 q 8 q 



+ -(-2 + 16z.-7z, 2 )^ + -(4 + 3^ 



2 1 



(n ■ p) 2 1 , 9n 1 , 

L*L+ (i-jv + J) (6.18) 



By looking at the structures in Eq. (|6.18|) we deduce that the most general form of G2PN is 



C 2 PN(g,p) = {q-p) 



« 2 p 4 + -(/? 2 p 2 + 7 2 (n-p) 2 ) + ^ 
q q\ 



(6.19) 
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Inserting the Ansatz (|6.19 ), and the lPN-results, in Eq. (|6.12 ), we get again more equations 
than unknowns: 



— a.2 



v 3 9 

_ + _ z/ 2 = . 



9 v 2 
4a 2 + /5 2 -3 7 2 + -z/ + — = 0. 

] + f3 2 - 5 2 - 2v+lis 2 = 0, 
4 o 

1 r 7 ^ n 

-I + ^ + I^-T = - 



a 2 - & + o - ^ 2 = 

o 

9 



3 72 — 3v — —v 







-i + 2/3 2 + 25 2 + 3 72 -^ 2 = 



(6.20) 



As it should (in view of the work of the previous sections) one finds that all the redundant 
equations can be satisfied. The final, unique solutions for the coefficients a 2 ,/5 2 ,7 2 and <5 2 
are: 



a 2 



72 



z/ + 3z/ 2 



8v + 3v 2 



62 



2v - 5z/ 2 



1 - Iv + v 2 



(6.21) 



Finally, we give the explicit form of the canonical transformation between the coordinates 
(q,p) and (q',p') at the 2PN level (see Eq. (|OD): 



q - q = — 

c z 
1 

+ -r 



v \ q v j „2 



1+ o --n?'p (g-p) 

2/ g 2 



^1 + 77 



z/\ g 4 (qr • p) 2 



+ - 5 - - 
4 V 2 



z/^ q l p 2 + 3 ^ ^ _ P l (q-p) 



q 



q 



+ \{l-7v + v 2 ) ^ + ^l-v) q *p* + ^l + v)p*p 2 ( q .p) 
4 g z 8 2 



(6.22) 



ft - ^ = -? 



v\ Pi v 9 
2/ g 2 F iP 



1 , ^ ft (g ■ g) 
2/ g 3 



+ 



1 

4 



(-1+31/) p,p 4 + ^(3 + llz/)^-^ (3 + 



+ ^(-2-18„ + , 2 ) *^ + £(10-i, 



g 2 4 

ft (9 ■ P) P 2 



2/ g 



1 2 3 . ^ Qi(q • p 



9 



(6.23) 



Note that the 1/ — ► limit of Eq. ( |6.22| ) gives g'* = (1 + l/(2c 2 g)) 2 g* which is (as it should) 
the relation between "Schwarzschild" (g') and "Isotropic" (g) quasi-cartesian coordinates in 
a Schwarzschild spacetime. [In this case, ADM = Isotropic]. As a check on Eqs. (|6.22|) , 
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(|6.23 ) we have verified that (at the 2PN level) q' x p' coincides with q x p. [They should 



coincide exactly, when solving exactly Eqs. ( |K.8| ) with any (spherically symmetric) generating 
function G(q,p).} Let us quote, for completeness, the partial derivatives of the generating 
function G = c~ 2 Gipn + c~ 4 G 2 pn, that must be used to solve by successive iterations the 
exact equations ( |6.8|) and determine q' and p' in terms of q and p: 

d Gl M<l,p)_ » pip2+ ( 1 + »\ *-(!+") *(«;P\ (6.24) 



dq i 1 V 2/ q V 2/ q 

5Gi PN (g,p) v i 2 . /, . v\ <t i, s f aOK \ 
---g'p + 1 + - vp l (q-p), (6.25) 



% 2^ V 2/ g 



5G 2 PN(g,p) i n . o \ 4,^/0 r \ Pi p 3 ,Pi(q-p) 
dqi =-u(l + 3u) PiP +-(2-5u) — + - V (S + 3 V )—^- 



q 5 4 g 2 8 g 3 
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VII. INCLUSION OF RADIATION REACTION EFFECTS AND TRANSITION 

BETWEEN INSPIRAL AND PLUNGE 

In the preceding sections we have limited our attention to the conservative (time- 
symmetric) part of the dynamics of a two-body system, i.e. the one defined, at the 2PN 
level, by neglecting A r ^ c in Eq. ( |1 . 1| ) . We expect that the separation of the dynamics in 
a conservative part plus a reactive part, makes sense also at higher PN orders (though it 
probably gets blurred at some high PN level). However, there exists, at present, no algo- 
rithm defining precisely this separation. Anyway we shall content ourselves here to working 
at the 2.5PN level where this separation is well-defined, as shown in Eq. (|1 . 1|) . When dealing 
with the relative motion we find it convenient to continue using an Hamiltonian formalism. 
Schafer f2"0|| , [I^j, fll~8f has shown how to treat radiation reaction effects within the ADM 



canonical formalism. His result (at the 2.5PN level) is that it is enough to use as Hamilto- 
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nian for the dynamics of two masses a time- dependent Hamiltonian obtained by adding to 
the conservative 2PN Hamiltonian H 2 p^(qi, q 2 ,Pi,p^) the following "reactive" Hamiltonian 



^eac(qi,q 2 ,Pl,P 2 ;t) = -^ TreaC (t) 



p{ A + pi pi _ 1 G m2 (?i - ?2)(si ^ 4) 



13 



(7.1) 



2mi 2m 2 2 |qr 1 — qr 2 | 

where 

> TTreac/,\ _ _^ G d 3 Qjjjt) . . 

^ W " 5 C 5 ^ > 

Qij denoting the quadrupole moment of the two-body system 



(7.3) 



Note that /ig Treac in Eq. should be treated as a given, time- dependent external field, 

considered as being independent of the canonical variables q a , p a . In other words, when 
writing the canonical equations of motion q = dH tot /dp, p = —dH tot /dq, one should 
consider only the explicit q — p dependence appearing in the square bracket on the R.H.S. of 
Eq. (|7.1|) . After differentiation with respect to q and p one can insert the explicit phase-space 
expression of the third time derivative of Qijit) (obtained, with sufficient precision, by using 
the Newtonian- level dynamics, i.e. by computing a repeated Poisson bracket of Qij(q,p) 
with H Newton (q,p)). 

Finally, we propose to graft radiation-reaction effects onto the non-perturbatively re- 
summed conservative dynamics defined by our effective-action approach in the following 
way. The total Hamiltonian for the relative motion Q, P in ADM coordinates is 

H tot (Q, P; t) = H™r vcd (Q, P) + H^ C (Q, P; t) , (7.4) 

where the "improved 2PN" Hamiltonian is that defined by solving Eq. (|4.25|) for £ real = H^, 
i.e. 



h ^q,p) l 1+2v ^ mQ ,nnQ,p)) _ ^ (75) 

on the R.H.S. of which one must transform, by the canonical transformation discussed in 



Sec. [VI|, the (exact) effective Hamiltonian defined by Eq. (|6.3|) . In the latter, we propose to 



use our current best estimates of the effective metric coefficients A(Q'), B(Q'), namely 
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B(Q') = A-\Q' 



(GMX 

WQ'J 



( GM\ 

1 - 6v W) 



(7.6) 



On the other hand the "reactive" contribution to the total Hamiltonian ( |7.4j ) is the center 
of mass reduction (p l = —p 2 = P,Q = q 1 — q 2 ) of Eq. ([7.1|) . 

In terms of reduced variables (q = Q/GM,p = P/ p) and of the non-relativistic reduced 
Hamiltonian, = (H^ al — Mc 2 )/ fi, our proposal reads 



H™(q,p; t) = if r ^ impr ° ved (g,p) + H- C (q,p; t) 



with 



NR improved / \ ^ 
real (?» PJ = ~ 



1 + 2z/ 



— H e s{q'{q,p),p > {q,p)) - 1 



where H e ff(q',p') is defined by inserting ( |7.6|) into Eq. ( |6.5|) , and with 



H^ c (q,p;t) = -^ c (t) 



1 ( ■ 1 gV 



(7.7) 



(7.8) 



(7.9) 



TTreac 



4 i/ 



-4(p l n 3 + p J n l ) + 6 n l n j (n ■ p) + - (n ■ p) 5 lj 

3 



(7.10) 



5 c 5 g 2 

where n % = <f/g. As explained above, the quantity hj^ reac {t) should not be differentiated 
with respect to q and p when writing the equations of motion 

_ = dHZ impmved (q,P) dH™*{q,p; ^ Trcac (t)) 
q dpi dpi 

q £NRimproved (?) p) Q j^reac^ p; h^ eac (t)) 



Pi = - ™* rl ■ ~^Ti — ■ (7-11) 

When inserting, after differentiation, Eq. ( |7.10| ), the equations of motion ( |7.11| ) become an 
explicit, autonomous (time-independent) evolution equation in phase space: x = f(x) where 
x = (q l ,Pi). From the study in Sec. above of the circular orbits defined by the exact, non- 
perturbative Hamiltonian H c s, we expect that the combined dynamics ( |7.11| ) will exhibit a 
transition from inspiral to plunge when q — \q\ (which decreases under radiation damping) 
reaches the image in the q — p phase space of the ISCO, studied above in q', p' coordinates. 
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FIG. 5. Inspiraling circular orbits in (q 1 ,p') coordinates including radiation reaction effects for 
v = 0.1 (left panel) and v = 1/4 (right panel). The location of the ISCO and of the horizon are 
indicated. 



We have in mind here quasi-circular, inspiraling orbits (circularized by radiation reaction), 
though, evidently, our approach can be used to study all possible orbits. We further expect 
that, when v <C 1 the inspiral will be very slow (the reaction Hamiltonian being proportional 
to v, see Eq. ( |7.10| )) and therefore the transition to plunge will be quite sharp, and well 
located at the ISCO. When v = 1/4 the radiation reaction effects are numerically smallish, 
but not parametrically small at the ISCO, and the transition to plunge cannot be expected 
to be very sharp. These expected behaviors are illustrated in Fig. [|. 

For simplicity, we have computed the orbits exhibited in these figures in q' space, neglect- 
ing the (formally 3.5PN) effect of the (q,p) —>■ [q',p') transformation on the reactive part of 
the equations of motion. [Thanks to the canonical invariance of the Hamilton equations of 
motion, the crucial conservative part of the evolution in q',p' space is simply obtained from 
the Hamiltonian H™ improved (q' , p') defined by keeping the variables q' and p' on the R.H.S. 
ofEq. Q.] 

Let us finally mention another possibility for incorporating radiation reaction effects 
directly in the effective one-body dynamics. In the q — p coordinates the (2.5PN) reaction 
Hamiltonian ( [7.1Q can be simply seen as due to perturbing the Euclidean metric = 5ij 
appearing in the lowest order Newtonian Hamiltonian (q l ab = q\ — q\) 

90 Pat Paj Gm a m b 
a<b 



^Newtonian ^a, Pa J ~ 2^ /_0 ~i J U/2 ' [ ' } 

a<b Wij lab lab) ' 
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by taking into account the near zone radiative field: 

9ij * 9% + hj^(t) , * g« - h^(t) . (7.13) 

By mapping back (through our (qp) «-> (q'p 1 ) link) the metric perturbation /i^ Treac onto 
the effective problem, one might try to incorporate reaction effects by defining a suitable 
"reactive" perturbation of our effective metric: 

M?') =<(<?') + * reac <(<?')■ (7-14) 

This approach might be useful for trying to go beyond the 2.5PN level discussed here and to 
define a "re-summed" version of reaction effects. Alternatively, if one has at one's disposal a 



more complete PN-expanded reactive force expressed in the original q coordinates ||32|| , one 
can, following the strategy proposed in Eq. ( |7.4| ), graft this improved (perturbative) reactive 
force onto the non-perturbatively improved conservative force defined by mapping back our 
effective dynamics onto the q coordinates. 

VIII. CONCLUSIONS 

We have introduced a novel approach to studying the late dynamical evolution of a 
coalescing binary system of compact objects. This approach is based on mapping (by a 
canonical transformation) the dynamics of the relative motion of a two-body system, with 
comparable masses mi, rri2, onto the dynamics of one particle of mass /i = mi 1712/(1711 +1712) 
moving in some effective metric ds e g. When neglecting radiation reaction, the mapping rules 
between the two problems are best interpreted in quantum terms (mapping between the 
discrete energy spectrum of bound states). They involve a physically natural transformation 
of the energy axis between the two problems, stating essentially that the effective energy of 
the effective particle is the energy of particle 1 in the rest-frame of particle 2 (or reciprocally), 
see Eq. (|4.26|) . The usefulness of this energy mapping was previously emphasized both in 
quantum two-body problems [[!]], and in classical ones [ ^jj . 



Starting from the currently most accurate knowledge of two-body dynamics |J, [0], we 
have shown that, when neglecting radiation reaction, our rules uniquely determine the effec- 
tive metric g e ^{q') in which the effective particle moves. This metric is a simple deformation 
of a Schwarzschild metric of mass M = mi+1712, with deformation parameter v = fi/M. Our 
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suggestion is then to define (as is done in quantum two-body problems [|TJ, 0) a particular 
non-perturbative re-summation of the usual, badly convergent, post-Newtonian-expanded 
dynamics by considering the dynamics defined by the effective metric as exact. This def- 
inition leads, in particular, to specific predictions for the characteristics of the innermost 
stable circular orbit (ISCO) for comparable-mass systems. In agreement with some previous 
predictions (notably one based on Pade approximants PHI), but in disagreement with the 
predictions of the "hybrid" approach of Ref. |[22| , we predict an ISCO which is more tightly 
bound than the usual test-mass-in-Schwarzschild one. The invariant physical characteristics 
of our predicted ISCO are given in Eqs. (|5.31|) and (|5.32|) , see also Tab. |. Note in par- 
ticular that the binding energy at the ISCO is robustly predicted to be ~ -1.5%Mc 2 
(for equal- mass systems; v = 1/4), while the orbital frequency at the ISCO is numerically 
predicted to be (again for v = 1/4) 

/ isco = 2372Hertz (Mq 

J V M 

Note that this corresponds to ~ 847 Hertz for (1.4M , 1.4M ) neutron star systems. 

We have argued, by studying the effects of higher (time-symmetric) post-Newtonian 
contributions, that our predictions for the characteristics of the ISCO are rather robust 
(especially when compared to the scatter of previous predictions). See Fig. |] and Tab. |. 
We note, however, that the knowledge of the 3PN dynamics (currently in progress [ 1191] , 
3q| ) would significantly reduce the present (2PN-based) uncertainty on the knowledge of 



the effective metric. 

The coordinate separation, in effective Schwarzschild coordinates, corresponding to the 
ISCO is Q' = R ~ 5.72GM/c 2 , i.e. ~ 23.6 km for a (1.4M , 1.4M ) neutron star system 
(from our canonical transformation ( |6.8D , this corresponds to an ADM-coordinate relative 
separation of Q ~ 4.79 GM/c 2 ). This value is near the sum of the nominal radii of (iso- 



lated) neutron stars for most nuclear equations of state |34| . This suggests that the inspiral 
phase of coalescing neutron star systems might terminate into tidal disruption (or at least 
tidally-dominated dynamics) without going through a well-defined plunge phase. Fully rel- 
ativistic 3D numerical simulations are needed to investigate this question. We note that a 
positive aspect of having (as predicted here) a rather low ISCO is that the end of the inspiral 
phase might well be very sensitive to the nuclear equation of state, so that LIGO/VIRGO 
observations might teach us something new about dense nuclear matter. 
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Finally, we have proposed two ways of adding radiation reaction effects to our effective 
one-body dynamics. The most straightforward one consists in directly combining radiation 
effects determined in the real two-body problem with the non-perturbative conservative dy- 
namics (which, in particular, features a dynamical instability at our ISCO) obtained by 
mapping the effective dynamics onto some standard (ADM or harmonic) two-body coordi- 
nate system: see Eq. (|7.7|) . A more subtle approach, which needs to be further developed, 
would consist in adding radiation reaction effects at the level of the effective metric itself, see 
Eq. ( 7.14 ). We have illustrated in Fig. |^ the transition from inspiral to plunge implied by (an 
approximation to) Eq. (|7. 7|) . In principle, this transition, and in particular the frequency at 
the ISCO, will be observable in gravitational wave observations of systems containing black 
holes. 

We hope that the approach presented here will also be of value for supplementing nu- 
merical relativity investigations. Indeed, our main (hopeful) claim is that the effective one- 
body dynamics is a "good" non-perturbative re-summation of the standard post-Newtonian- 
expanded results. Therefore, it gives a simple way of boosting up the accuracy of many 
PN-expanded results. [We leave to future work a more systematic analysis of the extension 
of our approach to higher post-Newtonian orders.] Effectively, this extends the validity of 
the post-Newtonian expansions in a new way (e.g., different from Pade approximants^). In 
particular, our results could be used to define initial conditions for two-body systems very 
near, or even at the ISCO, thereby cutting down significantly the numerical work needed to 
evolve fully relativistic 3D binary-system simulations. 

As a final remark, let us note that many extensions of the approach presented here are 
possible. In particular, the addition of the (classical) spin degrees of freedom to the effective 
one-body problem (in the effective metric and/or in the effective particle) suggests itself as 
an interesting issue (with possibly important physical consequences). 



7 It should be, however, possible to combine the effective one-body approach with Pade approxi- 
mants, thereby defining an even better scheme. 
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APPENDIX A: 



In this appendix we determine, at the 2PN level and in the Schwarzschild gauge, the 
effective metric 



ds 2 eS = -A(R) c 2 dt 2 + B(R) dR 2 + R 2 (d6 2 + sin 2 6 dip 2 ) 



(Al) 



A(R) = l + ^ + 4^+ " 3 



(A2) 



c 2 R ' c A R 2 ' c 6 i? 3 ' ~ ' c 2 R ' c 4 i? 2 ' 

when requiring simultaneously that: a) the energy levels of the "effective" and "real" 
problems coincide modulo an overall shift, i.e. £o{No,Jo) — ^reai(-A/", J) — Co, with 
c = M c 2 — m c 2 , Jo = J and Mo = M and b) the effective metric depends only on 
mi and m%. In this case, as anticipated in Sec. fVj, we will see that it not possible to satisfy 
the condition m = \i. 

The radial action I R (£o, Jo) of the "effective" description is 

2" 
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where f^ 111 = £ — m c 2 , a = GM m , 



(A3) 
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C = ai , 

4 64 



^ _ a 2 a 2 ai&i 



i? = — a 2 
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(A4) 



and we have introduced the dimensionless coefficients: 

^ a,- - h 



Cl; 



(GM ) 



(GM y ■ 



(A5) 
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We define the mass M used to scale the coefficients Oj and 6j by requiring % = —2 (i.e. 
a x = -2GM ). Identifying Eq. with the radial action /£>(£ NR , J) of the "real" problem, 
i.e. 
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NR 
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c 4 J 
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NR 



/iC z 



(A6) 



c 2 ^ 

where a = GM[i and £ NR = £ rea i— M c 2 , yields six equations to be satisfied. The requirement 
a) above implies the simple identification of the variables entering Eqs. ( |A3| ) and ( |A6| ) : 
£q R = £ NR , Jo = J, Ir — Ir- The explicit form of the equations stating that Am^ 2 a 
(OPN level), B tuq 1/2 a Q , D a 2 Q ( 1PN level ) and C m " 3/2 a , Ea 2 /m and Fa 4 (2PN level) 
in Eq. (|A3|) coincide with the analogous coefficients in Eq. ( |A6|) yields: 



-1 /2 




— v H i/' 

16 32 



-3i/ 



QT 




m^ 2 ao = /i 1 ^ 2 a , (A7) 

(A8) 

(A9) 
(AlO) 

(All) 
(A12) 

It is to be noted that if we impose = fi and GMq = GM (so that ao = a) we 
get an incompatibility at the 2PN level. Indeed, Eq. ([AT]) is satisfied and we can solve 
Eqs. (|A8|),([^9D in terms of the lPN-coefficients b\ and a2, but then the 2PN-equation 
( |A10| ), which contains only 6 1; is not satisfied. [This problem is due to the fact that 
we have more equations than unknowns.] Hence, we are obliged to relax the constraint 
m = /i. Let us introduce the parameter £, defined by m = yu£~ 2 . Eq. (|A7|) then gives 
GMq = GM C, 3 . Note that we are crucially using here the fact that the Newton-order energy 
levels £ NR = — m Q al/ (2Afo) + (D(c~ 2 ) do not depend separately on mo and ao = G M mo, 
but only on the combination m a 2 , = G 2 M5 2 m^. Solving the lPN-level Eqs. (|A8| ) , ( |A9|) we 
then get: 
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'i-TTil 15 - 7 ? 2 -"). S 2 = 775 (-9 + 9C 2 - v) , (A13) 



i_ (15 -7^-,), S 2 = jl 2 
while the 2PN-level Eq. ( |A10|) gives a quadratic equation in £ 2 which fixes uniquely its value 



(as well as that of the positive parameter £), namely 



m 5 



-15 + v + 2v / 2V50 + 15z/ + 2z/ 2 



(A14) 



Finally, the remaining 2PN equations ( All ) and ( |A12|) determine the coefficients of the 
effective metric at the 2PN level 

7 ' (ll85 -978£ 2 + 49£ 4 -414i/ + 14£ 2 z/ + z/ 2 ) , (A15) 



64 £ 2 
1 



a 3 = — — ■ (-289 + 402 £ 2 - 113 £ 4 + 158 v + 50£ 2 v - v 2 ) . (A16) 
64 £ 4 

The complexity of the results ( |A13| )-( |Al6| ), compared to the simplicity of our preferred 
solution ( |5.6| )-( |5T8"D , convinced us that the requirement a) above should be relaxed. Also, 
it seems fishy to have an effective mass tuq which differs from /i even in the non-relativistic 
limit c — > oo. Finally, it is not evident that this method can be generalized to higher 
post-Newtonian orders (where more redundant equations will have to be satisfied). 



APPENDIX B: 

In this appendix we describe an alternative, more formal method to map the "effective" 
one-body problem onto the "real" two-body one. We work in the Schwarzschild gauge. Here 
we require simultaneously that: a) the energy levels of the "effective" and "real" descriptions 
coincide modulo an overall shift, i.e. 8q(No, Jo) = S Tea \(J\f, J) — Cq, with Co = M c 2 — m c 2 , 
Jo = J and Ao = Af and b) the effective mass m is equal to the reduced mass \i = 

m2/(?7ii + m 2 ). Introducing the dimensionless quantities: 

rO rreal pNR pNR 
70 _ 1 R p-eal _ 1 R rp _ °0 771 _ '-'real 
J-R = , 1r = , &0 = , -^real = , 

ao a mo fi 

jo = —, j = -, (Bl) 

ao a 

where ao = GM mo and a = GM \x = Gm 1 m 2 , we can re- write the radial action for the 
"effective" problem, Eq. ( |3.13| ), in the form: 
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and where we have used, as above, the scaled metric coefficients 
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Identifying i^(-Eo> j ) with the analogous expression for the "real" problem, 
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and imposing E = £7 rea i, wiq = /i, «o = we get more equations to be satisfied than 
unknowns, 

— 01 = 1 , 
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Note that Eqs. (|B7|) and ( p9|) depend only on ai and b\, and cannot both be satisfied. To 
solve this incompatibility we consider here the possibility that the various coefficients that 
appear in the effective metric depend on the energy. Namely, at the 2PN level we consider 
the following expansions 
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h(E ) =b? + b? (^) , W)=^. (B15) 



The introduction of an energy dependence in the coefficients a^, hi reshuffles the c~ 2 expansion 
of Eq. (p2|) and modifies the equations (p6|)-( pTTD to be satisfied. It is easy to see that the 
flexibility introduced by the new coefficients af n \ bf^ allows one to solve in many ways 
the constraints to be satisfied. The simplest solution is obtained by requiring that the 
energy-dependence enters only in ai(E ) and only at the 2PN level: 

af > = , af = , bf ] = , (BI6) 



because in this case only Eq. (|B~^ ) gets modified. Indeed, it is straightforward to derive the 
new equation replacing (|B9|) : 



19^ ( o) &f a } 35 15 3 2 

a\ ' + — — = 1 v^ v 2 . (B17) 
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Hence, from Eqs. (B6)~(B£) we obtain the effective metric coefficients at the 1PN level: 



^ = -2, 3f> = -£, b? = \(8-,), (B18) 



while the 2PN-equations (|B17|) and flBTlD , (|BT2|) give: 



i? = -^(32 + 3u), 4 0) = ^ (208 - i/) , 6< 0) = ^(256 -400z/ + z/ 2 ). (B19) 



Again this solution is more complex than our preferred solution ( |5.6|) (|57q) . Moreover we 
think that the assumption of an energy dependence in the effective metric introduces a 
conceptual obscurity in the entire approach: Indeed, one should introduce two separate 
(effective) energies: the energy parameter E^ appearing explicitly in g^,, and the conserved 
energy Eq of some individual geodesic motion in g^l(E^). They can only be identified, a 
posteriori, for each specified geodesic motion. This makes it also quite difficult to incorporate 
radiation reaction effects. 
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Finally, one can require that the effective metric does not depend on the energy, but that 
the effective mass m depends on E . One then finds the solution 



m (E ) = /I 



(B20) 



with a corresponding effective metric defined by the energy-independent part a\°\ bf^ of the 
solution above. The objections of complexity and conceptual obscurity raised above also 
apply to this energy-dependent effective-mass solution. 
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